IN SILICO SCREENING OF THE POTENTIAL SARS-CoV-2 INHIBITORS BLOCKING THE HR1 TRIMER OF THE CORONAVIRUS PROTEIN S

A virtual library of biologically active molecules has been formed and in silico screening has been carried out for identification of small-molecule chemical compounds - potential SARS-CoV-2 inhibitors able to bind to the HR1 trimer of the protein S and to block the formation of a six-helix bundle 6-HB, which is critical for the virus-cell fusion and viral in-fectivity. Molecular modeling methods were used to evaluate the binding affinity of the identified compounds to the HR1 trimer of the protein S. As a result, 12 molecules exhibiting the high binding affinity to this functionally important region of the virus were found. The data obtained indicate the promise of using these compounds in the development of new antiviral drugs presenting SARS-CoV-2 fusion inhibitors that can block the virus entry into the host cell.

Novel Engineered SARS-CoV-2 HR1 Trimer Exhibits Improved Potency and Broad-Spectrum Activity against SARS-CoV-2 and Its Variants.

The ongoing pandemic of COVID-19, caused by SARS-CoV-2, has substantially increased the risk to global public health. Multiple vaccines and neutralizing antibodies (nAbs) have been authorized for preventing and treating SARS-CoV-2 infection. However, the emergence and spread of the viral variants may limit the effectiveness of these vaccines and antibodies. Fusion inhibitors targeting the HR1 domain of the viral S protein have been shown to broadly inhibit SARS-CoV-2 and its variants. In theory, peptide inhibitors targeting the HR2 domain of the S protein should also be able to inhibit viral infection. However, previously reported HR1-derived peptide inhibitors targeting the HR2 domain exhibit poor inhibitory activities. Here, we engineered a novel HR1 trimer (HR1MFd) by conjugating the trimerization motif foldon to the C terminus of the HR1-derived peptide. HR1MFd showed significantly improved inhibitory activity against SARS-CoV-2, SARS-CoV-2 variants of concern (VOCs), SARS-CoV, and MERS-CoV. Mechanistically, HR1MFd possesses markedly increased α-helicity, thermostability, higher HR2 domain binding affinity, and better inhibition of S protein-mediated cell-cell fusion compared to the HR1 peptide. Therefore, HR1MFd lays the foundation to develop HR1-based fusion inhibitors against SARS-CoV-2. IMPORTANCE Peptides derived from the SARS-CoV-2 HR1 region are generally poor inhibitors. Here, we constructed a trimeric peptide HR1MFd by fusing the trimerization motif foldon to the C terminus of the HR1 peptide. HR1MFd was highly effective in blocking transductions by SARS-CoV-2, SARS-CoV-2 variants, SARS-CoV, and MERS-CoV pseudoviruses. In comparison with HR1M, HR1MFd adopted a much higher helical conformation, better thermostability, increased affinity to the viral HR2 domain, and better inhibition of S protein-mediated cell-cell fusion. Overall, HR1MFd provides the information to develop effective HR1-derived peptides as fusion inhibitors against SARS-CoV-2 and its variants.

Structural conservation among variants of the SARS-CoV-2 spike postfusion bundle.

Variants of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) challenge currently available COVID-19 vaccines and monoclonal antibody therapies due to structural and dynamic changes of the viral spike glycoprotein (S). The heptad repeat 1 (HR1) and heptad repeat 2 (HR2) domains of S drive virus­host membrane fusion by assembly into a six-helix bundle, resulting in delivery of viral RNA into the host cell. We surveyed mutations of currently reported SARS-CoV-2 variants and selected eight mutations, including Q954H, N969K, and L981F from the Omicron variant, in the postfusion HR1HR2 bundle for functional and structural studies. We designed a molecular scaffold to determine cryogenic electron microscopy (cryo-EM) structures of HR1HR2 at 2.2­3.8 Å resolution by linking the trimeric N termini of four HR1 fragments to four trimeric C termini of the Dps4 dodecamer from Nostoc punctiforme. This molecular scaffold enables efficient sample preparation and structure determination of the HR1HR2 bundle and its mutants by single-particle cryo-EM. Our structure of the wild-type HR1HR2 bundle resolves uncertainties in previously determined structures. The mutant structures reveal side-chain positions of the mutations and their primarily local effects on the interactions between HR1 and HR2. These mutations do not alter the global architecture of the postfusion HR1HR2 bundle, suggesting that the interfaces between HR1 and HR2 are good targets for developing antiviral inhibitors that should be efficacious against all known variants of SARS-CoV-2 to date. We also note that this work paves the way for similar studies in more distantly related viruses.

Highlight of potential impact of new viral genotypes of SARS-CoV-2 on vaccines and anti-viral therapeutics.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causal agent of the coronavirus disease (COVID-19) pandemic, has infected millions of people globally. Genetic variation and selective pressures lead to the accumulation of single nucleotide polymorphism (SNP) within the viral genome that may affect virulence, transmission rate, viral recognition and the efficacy of prophylactic and interventional measures. To address these concerns at the genomic level, we assessed the phylogeny and SNPs of the SARS-CoV-2 mutant population collected to date in Iran in relation to globally reported variants. Phylogenetic analysis of mutant strains revealed the occurrence of the variants known as B.1.1.7 (Alpha), B.1.525 (Eta), and B.1.617 (Delta) that appear to have delineated independently in Iran. SNP analysis of the Iranian sequences revealed that the mutations were predominantly positioned within the S protein-coding region, with most SNPs localizing to the S1 subunit. Seventeen S1-localizing SNPs occurred in the RNA binding domain that interacts with ACE2 of the host cell. Importantly, many of these SNPs are predicted to influence the binding of antibodies and anti-viral therapeutics, indicating that the adaptive host response appears to be imposing a selective pressure that is driving the evolution of the virus in this closed population through enhancing virulence. The SNPs detected within these mutant cohorts are addressed with respect to current prophylactic measures and therapeutic interventions.

The fatty acid site is coupled to functional motifs in the SARS-CoV-2 spike protein and modulates spike allosteric behaviour.

The SARS-CoV-2 spike protein is the first contact point between the SARS-CoV-2 virus and host cells and mediates membrane fusion. Recently, a fatty acid binding site was identified in the spike (Toelzer et al. Science 2020). The presence of linoleic acid at this site modulates binding of the spike to the human ACE2 receptor, stabilizing a locked conformation of the protein. Here, dynamical-nonequilibrium molecular dynamics simulations reveal that this fatty acid site is coupled to functionally relevant regions of the spike, some of them far from the fatty acid binding pocket. Removal of a ligand from the fatty acid binding site significantly affects the dynamics of distant, functionally important regions of the spike, including the receptor-binding motif, furin cleavage site and fusion-peptide-adjacent regions. Simulations of the D614G mutant show differences in behaviour between these clinical variants of the spike: the D614G mutant shows a significantly different conformational response for some structural motifs relevant for binding and fusion. The simulations identify structural networks through which changes at the fatty acid binding site are transmitted within the protein. These communication networks significantly involve positions that are prone to mutation, indicating that observed genetic variation in the spike may alter its response to linoleate binding and associated allosteric communication.

Analysis of SARS-COV2 spike protein variants among Iraqi isolates.

The ongoing pandemic of COVID-19 caused by the SARS-COV2 virus has triggered millions of deaths around the globe. Emerging several variants of the virus with increased transmissibility, the severity of disease, and the ability of the virus to escape from the immune system has a cause for concerns. Here, we compared the spike protein sequence of 91 human SARS CoV2 strains of Iraq to the first reported sequence of SARS-CoV2 isolate from Wuhan Hu-1/China. The strains were isolated between June 2020 and March 2021. Twenty-two distinct mutations were identified within the spike protein regions which were: L5F, L18F, T19R, S151T, G181A, A222V, A348S, L452 (Q or M), T478K, N501Y, A520S, A522V, A570D, S605A, D614G, Q675H, N679K, P681H, T716I, S982A, A1020S, D1118H. The most frequently mutations occurred at the D614G (87/91), followed by S982A (50/91), and A570D (48/91), respectively. In addition, a distinct shift was observed in the type of SARS-COV2 variants present in 2020 compared to 2021 isolates. In 2020, B.1.428.1 lineage was appeared to be a dominant variant (85%). However, the diversity of the variants increased in 2021, and the majority (73%) of the isolated were appeared to belong to B.1.1.7 lineage (VOC/alpha variants). To our knowledge, this is the first major genome analysis of SARS-CoV2 in Iraq. The data from this research could provide insights into SARS-CoV2 evolution, and can be potentially used to recognize the effective vaccine against the disease.

COVID-19: the CaMKII-like system of S protein drives membrane fusion and induces syncytial multinucleated giant cells.

The SARS-CoV-2 S protein on the membrane of infected cells can promote receptor-dependent syncytia formation, relating to extensive tissue damage and lymphocyte elimination. In this case, it is challenging to obtain neutralizing antibodies and prevent them through antibodies effectively. Considering that, in the current study, structural domain search methods are adopted to analyze the SARS-CoV-2 S protein to find the fusion mechanism. The results show that after the EF-hand domain of S protein bound to calcium ions, S2 protein had CaMKII protein activities. Besides, the CaMKII_AD domain of S2 changed S2 conformation, facilitating the formation of HR1-HR2 six-helix bundles. Apart from that, the Ca2+-ATPase of S2 pumped calcium ions from the virus cytoplasm to help membrane fusion, while motor structures of S drove the CaATP_NAI and CaMKII_AD domains to extend to the outside and combined the viral membrane and the cell membrane, thus forming a calcium bridge. Furthermore, the phospholipid-flipping-ATPase released water, triggering lipid mixing and fusion and generating fusion pores. Then, motor structures promoted fusion pore extension, followed by the cytoplasmic contents of the virus being discharged into the cell cytoplasm. After that, the membrane of the virus slid onto the cell membrane along the flowing membrane on the gap of the three CaATP_NAI. At last, the HR1-HR2 hexamer would fall into the cytoplasm or stay on the cell membrane. Therefore, the CaMKII_like system of S protein facilitated membrane fusion for further inducing syncytial multinucleated giant cells.

Structural stability predictions and molecular dynamics simulations of RBD and HR1 mutations associated with SARS-CoV-2 spike glycoprotein.

The COVID-19 pandemic is caused by human transmission and infection of Severe Acute Respiratory Syndrome Corona Virus-2 (SARS-CoV-2). There is no trusted drug against the virus; hence, efforts are on discovering novel inhibitors for the virus. The entry of a SARS-CoV-2 virus particle into a host cell is initiated by its spike glycoprotein and host Angiotensin-Converting Enzyme 2 (ACE2) receptor interaction. Spike glycoprotein domains, namely, the Receptor Binding Domain (RBD) and Heptad Repeat (HR) domains, are essential for this activity. We have studied the impact of mutations such as A348T, N354D, D364Y, G476S, V483A, S494D in the RBD (319-591), and S939F, S940T, T941A, S943P (912-984) in the HR1 domains of spike glycoprotein. Summarily, we utilized the computational screening algorithms to rank the deleterious, damaging and disease-associated spike glycoprotein mutations. Subsequently, to understand the changes in conformation, flexibility and function of the spike glycoprotein mutants, Molecular Dynamics (MD) simulations were performed. The computational predictions and analysis of the MD trajectories suggest that the RBD and HR1 mutations induce significant phenotypic effects on the pre-binding spike glycoprotein structure, which are presumably consequential to its binding to the receptor and provides lead to design inhibitors against the binding.Communicated by Ramaswamy H. Sarma.

Unlocking COVID therapeutic targets: A structure-based rationale against SARS-CoV-2, SARS-CoV and MERS-CoV Spike.

There are no approved target therapeutics against SARS-CoV-2 or other beta-CoVs. The beta-CoV Spike protein is a promising target considering the critical role in viral infection and pathogenesis and its surface exposed features. We performed a structure-based strategy targeting highly conserved druggable regions resulting from a comprehensive large-scale sequence analysis and structural characterization of Spike domains across SARSr- and MERSr-CoVs. We have disclosed 28 main consensus druggable pockets within the Spike. The RBD and SD1 (S1 subunit); and the CR, HR1 and CH (S2 subunit) represent the most promising conserved druggable regions. Additionally, we have identified 181 new potential hot spot residues for the hSARSr-CoVs and 72 new hot spot residues for the SARSr- and MERSr-CoVs, which have not been described before in the literature. These sites/residues exhibit advantageous structural features for targeted molecular and pharmacological modulation. This study establishes the Spike as a promising anti-CoV target using an approach with a potential higher resilience to resistance development and directed to a broad spectrum of Beta-CoVs, including the new SARS-CoV-2 responsible for COVID-19. This research also provides a structure-based rationale for the design and discovery of chemical inhibitors, antibodies or other therapeutic modalities successfully targeting the Beta-CoV Spike protein.

Minireview of progress in the structural study of SARS-CoV-2 proteins.

A severe form of pneumonia, named coronavirus disease 2019 (COVID-19) by the World Health Organization, broke out in China and rapidly developed into a global pandemic, with millions of cases and hundreds of thousands of deaths reported globally. The novel coronavirus, which was designated as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), was identified as the etiological agent of COVID-19. On the basis of experience accumulated following previous SARS-CoV and MERS-CoV outbreaks and research, a series of studies have been conducted rapidly, and major progress has been achieved with regard to the understanding of the phylogeny and genomic organization of SARS-CoV-2 in addition its molecular mechanisms of infection and replication. In the present review, we summarized crucial developments in the elucidation of the structure and function of key SARS-CoV-2 proteins, especially the main protease, RNA-dependent RNA polymerase, spike glycoprotein, and nucleocapsid protein. Results of studies on their associated inhibitors and drugs have also been highlighted.